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ABSTRACT We investigated the structure of the hydrophobic domain of the severe acute respiratory syndrome E protein in 
model lipid membranes by x-ray reflectivity and x-ray scattering. In particular, we used x-ray reflectivity to study the location of an 
iodine-labeled residue within the lipid bilayer. The label imposes spatial constraints on the protein topology. Experimental data 
taken as a function of protein/lipid ratio P/L and different swelling states support the hairpin conformation of severe acute respiratory 
syndrome E protein reported previously. Changes in the bilayer thickness and acyl-chain ordering are presented as a function 


of P/L, and discussed in view of different structural models. 


INTRODUCTION 


The severe acute respiratory syndrome (SARS) that broke 
out in April 2003 is a newly identified infectious disease 
(1-5). The coronavirus (SARS-CoV) has been identified as 
the primary causative agent for SARS. Sequence analysis 
reveals the phylogeny of SARS-CoV, showing characteristic 
features of a coronavirus. At the same time it belongs to a 
new group that is sufficiently different from known corona- 
viruses (6,7). Coronaviruses have four important viral genes 
with different structural proteins: a spike glycoprotein (S), a 
small envelope protein (E), a matrix glycoprotein (M), and a 
nucleocapsid protein (N). In this article, we address the 
structure of the transmembrane domain of the SARS-CoV E 
protein in a model phospholipid membrane by x-ray scat- 
tering. 

S, M, and N proteins of different coronaviruses have been 
broadly studied for their important roles in receptor binding 
and virion budding. The significance of the E protein has 
been realized only much more recently (8). This membrane- 
bound constituent of the virion was not immediately rec- 
ognized as a viral structural protein, owing to its small size 
(~10 kDa) and its very low abundance relative to the M, N, 
and S proteins. Today, coronavirus E proteins are known to 
play an important role in viral morphogenesis. Coexpression 
of mouse hepatitis virus (MHV) E and M proteins results in 
the production of virus-like particles (9,10), indicating that 
neither the nucleocapsid nor the viral spike are needed for 
viral budding. The interaction between the two proteins 
relevant for the budding process is thought to take place in 
pre-Golgi compartments whereby the cytoplasmic domains 
of the two proteins interact (11). During the expression of E 
protein the Golgi apparatus changes its morphology dramat- 
ically (12), explaining in part E protein’s ability to induce 
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apoptosis (13,14). Expression of M protein from several 
coronaviruses on its own did not produce virus-like particles 
(9,10,15-17). On the contrary, expression of MHV E protein 
on its own caused the release of vesicles containing E pro- 
tein, highlighting the important role of this small protein. 
This pattern seems to be a general result for coronaviridea 
(15,16). 

E proteins are well conserved within each of the different 
groups of coronaviruses (18). Regarding the sequence, it 
is possible to make the following generalization for all E 
proteins: E proteins are all small proteins (~75 residues) 
with an unusually long hydrophobic stretch (25-30 amino 
acids), placed in between a hydrophilic N- and C-terminus, 
~8 and ~40 residues long, respectively. Note that the length 
of the hydrophobic segment of E proteins is significantly 
larger than the average length of a transmembrane a-helix, 
which is only ~21 residues (19). The SARS-CoV E (SARS- 
CoV E) protein is a 76-residue (NH;’-MYSFVSEETG- 
TLIVNSVLLFLAFVVFLL VTLAILTALRLCAYCCNI- 
VNVSLVKPTVYVYSRVKNLNSSEGVPDLLV-COO_ ) 
polypeptide. The region in bold type in the sequence represents 
the hydrophobic transmembrane domain (TMD) of the protein. 

In this article, we investigate the structure of dimyristoyl- 
phosphatidylcholine (DMPC) bilayers with reconstituted 
TMD of SARS-CoV E by x-ray reflectivity and grazing 
incidence scattering, as a function of peptide/lipid molar 
ratio P/L and the hydration level of the headgroups. Peptides 
were synthesized that encompassed the entire hydrophobic 
region of the protein (Glu-7 to Arg-38), as depicted in Fig. 
1 a. Ina previous study, we focused on the secondary struc- 
ture of TMD using Fourier transform infrared (FTIR) spec- 
troscopy, x-ray reflectivity, and molecular modeling (20). 
FTIR dichroism on oriented bilayers showed a highly helical 
backbone structure oriented perpendicular to the bilayer. At 
the same time, x-ray reflectivity analysis comparing an 
unlabeled peptide and a peptide labeled with iodine at po- 
sition Phe-23 in the center of the hydrophobic stretch 
allowed us to pinpoint the location of the central Phe-23 
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(a) Ribbon diagram of the structural model derived from the TMD of SARS-CoV E protein taken from Arbely et al. (20). (b) The sequence 


diagram of the SARS-CoV E protein around a pseudo point of symmetry (i-e., Phe-23) as two short helices forming the hairpin model. (c) Sketch illustrating 
a bilayer with incoming and reflecting beams in a reflectivity setup. a; and a¢ denote the angles between the sample surface and the incoming and reflected 
beams, respectively. g, denotes the momentum transfer in the z-direction, and qj the lateral momentum transfer. Also shown is a schematic of the SARS-CoV E 
protein embedded in the membrane bilayer. The dot represents the iodinated Phe-23 group. 


group in the bilayer electron-density profile. It was found to 
be located at a position displaced ~ 16-17 A from the bilayer 
center in the headgroup region. Here we give a more com- 
plete account of these experiments, complemented by addi- 
tional results obtained at a wider range of P/L, and also for 
three different levels of hydration. Furthermore, we inves- 
tigated changes in the acyl-chain ordering at high protein con- 
centrations. Finally, anomalous x-ray reflectivity was used to 
verify the location of the iodine label. 

Concluding from previous FTIR and x-ray results (20) for 
high P/L samples, a short hairpin conformation inserted 
perpendicular to the plane of the membrane was postulated 
for the sequence (see Fig. 1 5). The hairpin forms an in- 
version about a pseudocenter of symmetry (see Fig. 1 5) . 
Molecular modeling supported the model and showed that 
the two helices are likely to be stabilized by specific bonds, 
namely a salt bridge between Glu-8 and Arg-38, and a hy- 
drogen bond between a single asparagine amino acid in the 
two helices and the backbone amide group (see Fig. | a). An 
alternative model had been proposed by Shen and coworkers 
for the SARS-CoV E protein (21), based on CD spectra taken 
in aqueous solution. According to this model, the SARS- 
CoV E protein sequence forms a single transmembrane (TM) 
helix, and a short B-sheet segment forming a hydrogen bond 
with the lipid bilayer. This model cannot be reconciled with a 
Phe-23 position displaced from the bilayer center. It would 
also lead to a significant hydrophobic mismatch between 
hydrophobic amino acids and acyl chain thickness, since the 
average number of amino acids in a transmembrane helix 
needed to span the entire bilayer is only ~21. Furthermore, 
recent antibody binding results with an epitope-tagged MHV 
E protein are indicative of the protein traversing the lipid 


bilayer twice, whereby both termini of the protein reside in 
the virus lumen (22). This can well be understood, if the 
hairpin conformation is the general motif of E proteins. On 
the other hand, molecular dynamics simulation has been 
used recently as a test for evolutionary conservation using 
several coronavirus protein homologous sequences, and 
points to a transmembrane oligomer topology (23). Note that 
the x-ray study alone is not able to prove or disprove the 
hairpin conformation. However, it can provide important stru- 
ctural constraints, as well as the electron-density profile of the 
lipid bilayer as a function of P/L (see Fig. 1 c). 

This introduction is followed by the sections Materials and 
Methods, X-ray Reflectivity and Electron-Density Profiles, 
Anomalous Reflectivity, and Acyl Chain Ordering Induced 
by E Protein. Finally, the article closes with Summary and 
Conclusions. 


MATERIALS AND METHODS 
Materials 


Dimyristoyl-sn-glycero-3-phosphatidylcholine was purchased from Avanti 
Polar lipids (Alabaster, AL). The purity of DMPC is claimed to be 99%. 
Therefore the lipid was used without further purifications. Chloroform (Chl) 
and 1,1,1,3,3,3-hexafluor-2-propanol (HFI) (purity 99.8%) were purchased 
from Sigma (Schnelldorf, Germany). 


Peptide synthesis 


The peptide was synthesized and purified (residues 7-38) by standard solid- 
phase N-(9-fluorenyl) methoxycarbonyl (Fmoc) chemistry as described in 
Arbely et al. (20). Two different synthetic peptides were made: an unlabeled 
peptide and one that contains iodine at position 23 (i.e., phenylalanine) of the 
sequence. 


Biophysical Journal 90(6) 2038-2050 


2040 


Preparation of lipid-protein multilamellar stacks 


Lipids were used as purchased to prepare multilamellar bilayers of DMPC/ 
protein complexes following the procedure described by Seul and Sammon 
(24). The lipids were first dissolved in a solution of Chl/HFI (1:1, v/v) at a 
concentration of 20 mg/ml, whereas the proteins were dissolved in a solution 
of Chl/HFI (40:60%, v/v) at 2 mg/ml since the protein is more soluble at 
higher HFI concentrations. Varied amounts of the protein stock solution 
were then mixed with the DMPC stock solution at a final concentration of 
5 mg/ml, which yields the desired P/L ratio. Pure solvents were added to 
yield an identical final lipid concentration. The protein/lipid ratio P/L ranged 
from 1/500 to 1/7.5. The mixed solutions were spread on Silicon substrates, 
cleaned by two 15-min cycles of ultrasonic bath in methanol, followed by 
two 15-min cycles in ultrapure water (18 MO. cm, Millipore, Bedford, MA), 
and finally dried under a nitrogen stream. A droplet of 200 wl was then 
spread on the Si-wafer of typically 15 X 25 mm? positioned in an exactly 
horizontal plane. The spread solution was allowed to dry very slowly to 
prevent film rupture and dewetting. The samples were then exposed to high 
vacuum for 12 h to remove completely all the solvent traces. Afterward, the 
samples were rehydrated, yielding film thicknesses in the range of 
D ~2— 5pm. Such a procedure produces multilamellar stacks well aligned 
with respect to the substrate, with a typical mosaicity (orientational 
distribution) less than the instrumental resolution (i.e., 0.01°) (25). A very 
low mosaicity is a prerequisite in applying interface-sensitive x-ray 
scattering techniques for structural studies of solid-supported bilayers. To 
examine the sample quality, we used light microscopy in bright-field 
contrast (Olympus, Melville, NY; objective, Neofluar 10X/0.3) to image the 
samples after the deposition at different P/L. The samples were kept in the 
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fluid state at the same temperature as for the x-ray experiments. A significant 
effect of the concentration of the SCoV E protein on the morphology of the 
supported membranes was observed. Images recorded in the Ly state at 
temperature T = 45° and relative humidity R.H. = 98% are shown in Fig. 2. 
The samples were kept in a sealed temperature-controlled chamber (Julabo, 
Seelbach, Germany) with a water reservoir at the bottom for hydration. Note 
that multilamellar films are known to exhibit a pronounced domain struc- 
ture with a large variation in local film thickness. At the same time, the 
orientation of the bilayers is almost perfect, despite the limited lateral 
extension of the domains, leading to a patch-like morphology (see Fig. 2 a 
for the case of pure DMPC). 

A significant change in the domain size and texture is observed already at 
a relatively small concentration of SARS-CoV E protein. At P/L = 0.002, 
the formation of small irregularly shaped domains, which comprise a sig- 
nificant fraction of the bilayer surface area, are observed instead of the 
relatively large domain structures in pure DMPC films (see Fig. 2 b). This 
defect structure is interconnected, and, as the concentration is increased to 
P/L = 0.01, changes to a pattern with even smaller domains. The typical 
length scales of the defect structures remain small for P/L = 0.01 and 0.02 
(Fig. 2, c and d), until at higher concentrations (i.e., 0.05 and 0.10) it 
increases again. The pattern at P/L = 0.05 exhibits many individual 
structures that are not connected (Fig. 2 e). Finally, a star-like morphology 
with relatively large smooth areas evolves on the top of the lipid film (Fig. 
2 f). These changes are very reproducible for different samples and upon 
translation of the illuminated spot on the sample. The results show that in the 
La phase, the SARS-CoV E protein drastically affects the multilayer 
morphology, possibly by changing the line tension between the domains. 


FIGURE 2 L,-phase (i.e., T = 45°C and R.H. = 98%) 
protein/lipid film images by bright-field contrast for dif- 
ferent protein concentrations P/L. The defect structures and 
domain sizes change with P/L. Scale bar, 100 um. 
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X-ray reflectivity 


Before x-ray reflectivity measurements, the resulting multilamellar stacks 
were inserted in a closed temperature- and humidity-controlled chamber. 
The chamber consists of two concentric stainless steel cylinders, with kapton 
windows. The chamber temperature was maintained by a flow of oil 
connected to a temperature-controlled reservoir (Julabo). The temperature 
was measured close to the sample holder using a Pt100 sensor with thermal 
stability in the range of 0.02 K over several hours (26). The average 
temperature of the samples was kept at T = 45°C, well above that of the 
chain-melting transition. The samples were mounted in the inner cylinder of 
the chamber facing a humid atmosphere controlled by adding a salt to a 
water reservoir placed at the bottom of the cylinder (27). By changing the 
salt type one can vary the relative humidity from 11% to 100%. We used 
three types of salt, namely LiCl, NaCl, and K2SOy,, leading to R.H. = 11%, 
75%, and 98%, respectively. When using K2SO,, DMPC bilayers were 
typically swollen up to a repeat distance of d ~ 50 A in the L,-phase. At 
both R.H. = 75% and R.H. = 98%, the membranes were in the fluid L, state, 
whereas the R.H. = 11% curves were indicative of a gel phase. 

The reflectivity experiments were carried out on the bending magnet 
beamline D4 of the DORIS storage ring at the synchrotron radiation 
laboratory HASYLAB/DESY (Hamburg, Germany) using a photon energy 
of 11 keV (i.e., A = 1.13 A), set by a Si(111) monochromator. The chamber 
was mounted on the z-axis diffractometer with the samples oriented 
vertically. The reflectivity curves were measured with a fast scintillation 
counter (Cyberstar, Oxford Instruments, Eynsham, U.K.) using motorized 
collimating slits on both incident and reflected beam paths. The reflectivity 
curves were corrected for ring current, sample illumination, and diffuse back- 
ground (offset-scan). 

In the following text, we briefly repeat the principles of x-ray reflectivity 
and the Fourier synthesis (FS) method as tools to determine the electron 
density of the protein-lipid system (28). To record a reflectivity curve, the 
incident beam with wave vector k; has to be collimated to less than a few 
hundredths of a degree and directed on the sample at glancing incidence 
angle a;. The reflected intensity is then measured as a function of a; under 
specular conditions (e.g., at an exit angle a¢ = a@;), with the wave vector of 
the exit beam denoted by k¢. Thus, the momentum transfer of the elastic 
scattering g = kp — k, is always along q,, with the z axis parallel to the sample 
normal (Fig. 1). Typically, the reflectivity can be recorded over seven to 
eight orders of magnitude (after correction for diffuse scattering and back- 
ground), as measured in a so-called offset scan. To this end, the x-ray reflectivity 
in the semikinematic approximation from an interface characterized by electron- 
density profile p(z) between two media of electron densities p; and p2 is given 
by Braslau et al. (29) 
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= Ry (4) : ane) 


; 2 
22 | 
Ap, Oz liar 
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where Rx is the Fresnel reflectivity of the ideal (sharp) interface between the 
two media, q, is the scattering vector, and Ap; is the density contrast of the 
two mediums. In this formalism the interface normal is along the z axis and 
p(z) is the laterally averaged electron density. For the moment we ignore the 
effect of absorption, which can be accounted for by introducing an 
imaginary component of the wave vector. Then, the function Rp(g,) can be 
written in terms of the critical momentum transfer q. as lq, — qil/lqz + gil 
with g/? = Iq," — q,’l. The critical momentum transfer or the critical angle is 
directly related to the density contrast between the two media by 
de = 477/Asin(a,) & 4,/ aro Apyy, with rg denoting the classical electron 
radius. In this case, medium 1| corresponds to air and medium 2 corresponds 
to the solid substrate (i.e., silicon). The multilamellar stack of bilayers is 
modeled by an interface with a partially oscillatory density profile as 
described in Salditt et al. (28). Using the linearity of the integrand in Eq. 1, 
one can decompose it into parts, the first one accounting for the density 
increment at the substrate that does not depend on the multilamellar bilayers, 
and a second one that contains entirely the information about the bilayer 


R(qz) = Re(qz)|P (4) 
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stack. The second term can be broken up into two parts: the form factor f(q,) 
and the structure factor s(q,). The structure factor contains the parameters of 
the multilamellar stack. For an ideal one-dimensional stacking, s(g,) would 
be simply given by 


No 7 
=a eo (2) 
n=1 


where d is the periodicity and N is the total number of layers in the stack. 
The form factor, which characterizes the electron-density distribution, is 
defined in this context as 


*” Ap(z) 
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The electron-density profile of the bilayer can be defined in various ways. 
A practical parameterization is in terms of its Fourier coefficients, since the 
number of parameters can easily be adapted to the resolution of a particular 
experiment. The deviations from the average bilayer electron density po in 
terms of the first No Fourier coefficients f, can thus be written as 


pte) = 3 frcos(7™), a) 


n=1 


where the phases v, = +1 are reduced to positive/negative signs due to the 
mirror plane symmetry of the bilayer. In general, an educated guess of the 
phasing can be deduced from the basic bilayer profile or the data as dis- 
cussed in Salditt et al. (28). The coefficients f,, can be related to the integrated 
intensities under the reflectivity curve after application of (Lorentz) 
correction factors, e.g., g, | or gz ~. For highly oriented films, a correction 
factor of g, ' has been proposed in Tristram-Nagle et al. (30) and 
Gandhavadi et al. (31). However, in the absence of a rigorous derivation it 
remains unclear in which approximation a simple factor suffices to correct 
the raw data. It is clear that effects of mosaicity, absorption, and Fresnel 
reflectivity terms all influence the intensities of the Bragg peaks, leaving 
aside for the moment additional effects of the structure factor, e.g., due to 
thermal fluctuations. In Li et al. (32), an experimental approach is used to 
address these questions, where the profiles derived from full g, range fitting 
are compared to those computed by the FS method, with different correction 
factors. Here we use f,%./nl,, which can be regarded as an approximation. 
Using a more rigorous reflectivity analysis it was found that the approx- 
imation is quite good, in particular for the bilayer thickness, which is of 
importance here. Furthermore, a comparison of profiles for labeled and 
unlabeled proteins should be meaningful even if systematic errors persist in 
the profile resulting from the approximation in this data treatment. Finally, 
the electron density of the bilayer lipid membrane on an arbitrary scale (no 
absolute units) has been calculated form the measured peak intensities as 


p(z) = y V/nl,V_ COS (72). (5) 
n=1 


RESULTS AND DISCUSSION 
X-ray reflectivity and electron-density profiles 


Characteristic reflectivity curves of the noniodinated se- 
quence P/L series are shown in Fig. 3, shifted vertically for 
clarity, for three different swelling states (R.H.). Similar 
reflectivity curves were obtained for the analogous series 
with the iodinated sequence. The reflectivity is plotted as a 
function of the vertical momentum transfer q, after subtrac- 
tion of the diffuse scattering (offset scan), and after 
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FIGURE 3 
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Reflectivity curves of multilamellar DMPC/SARS-CoV E protein membranes at three distinct R.H. and constant T = 45°C. Different P/L ratios 


are presented for the noniodinated (upper panels) and iodinated (lower panels) SARS-CoV E protein. The curves are shifted vertically for clarity. The curves 
exhibit the typical pattern of a lamellar structure with well defined periodicities d and large vertical domain sizes as evidenced from the sharp peaks. The smaller 
peak heights for P/L = 1/10 reflect the increase of fluctuations and/or static disorder with P/L. All data has been corrected for illumination, diffuse background, 


and beam monitor. 


illumination correction. The curves show the typical features 
of highly oriented multilamellar films: the plateau of total 
reflection at small g,, and a set of sharp and intense 
equidistant Bragg peaks. The intensity and number of Bragg 
peaks decreases with P/L, indicative of peptide-induced 
lamellar disorder. The complete P/L series was measured at 
constant temperature T = 45°C, but at different swelling 
states, controlled by hydration from saturated salt solutions 
(see Materials and Methods), corresponding to nominal 
relative humidities of 98%, 75%, and 11%. The first two 
R.H. values correspond to the L, phase where the Bragg- 
peak intensities decay more homogenously. 

The gel-phase curves show a particularly interesting 
pattern upon addition of the protein, where odd reflection 
orders with the exception of n = | are significantly sup- 
pressed, but since we are more interested in the fluid state, we 
leave these curves aside for the moment. In the fluid state, the 
reflectivity curves depend in a systematic way on the peptide 
concentration. This dependency and the corresponding 
d-spacings are identical for the iodinated and noniodinated 
series, indicating that the presence of the label does not 
perturb the system. 

To evaluate the data, we used the Fourier synthesis method, 
using solely the integrated peak intensities to compute the 
density profile, rather than full g, range fits of the curves 
(28), as discussed in the previous section. The latter has a 
much larger potential for structural analysis and yields 
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density profiles in absolute units, but is also very time- 
consuming and more difficult to achieve. Although we are 
working toward this goal, we are not yet in possession of a 
comprehensive reflectivity fitting model and software for 
these films. A suitable model function must take into account 
effects of absorption, thermal fluctuations, static defects, and 
instrumental resolution, and yet keep the number of param- 
eters manageable. Therefore we turn here to the Fourier 
synthesis method, which is commonly used in the literature 
but which poses some serious problems, as discussed in Li 
(33). In particular, Fresnel-type reflectivity curves are not 
included in this simplified approach. Moreover, changes in 
the structure factor with P/L are falsely attributed to the form 
factor, since the structure factor is tacitly assumed to be that 
of the ideal lattice. More than seven lamellar reflectivity 
reflections were observed for the peptide-free bilayers, whereas 
only five or four orders persist in the presence of the peptide 
sequence, independent of iodination. This phenomenon is 
typical for many membrane active peptides or proteins and 
leads to a smoothing of the deduced bilayer profile. Of 
course, the local profile is not necessarily flatter for high P/L. 
Instead, this effect probably results from the increased 
lamellar disorder, since the determined profiles have to be 
regarded not as intrinsic profiles, but rather as convolutions 
of the intrinsic profiles with the distribution function of the 
bilayer position, which broadens with increasing lamellar 
disorder. Five orders of reflectivity data are clearly enough to 
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calculate the electron distribution of the peptide/bilayer with a 
sufficiently high resolution to determine structural quantities, 
as the distance between the headgroups, e.g., the distance 
between the two maxima corresponding to the phosphorus 
atoms d,,. A summary of the integrated peak intensities for 
98%, 75%, and 11% R.H. are presented in Table 1. 

The centrosymmetric electron-density profiles of the bi- 
layers containing noniodinated and iodinated SARS-CoV E 
protein obtained at R.H. = 98% for different P/L are shown 
in Fig. 4, on an arbitrary scale. The profiles were calculated 
using Eq. 6, with appropriate choice of phases (—, —, +, —, 
+, —, —). The curves have been normalized such that the 
area under the first Bragg peak is set to be 1, whereas the 
higher-order Bragg peaks were properly scaled with respect 
to the first peak intensity. Similar curves were obtained at 
75% and 11% R.H. The well known interpretation of the 
profiles is as follows: the two mean peaks of p(z) on either 
side of the figure correspond to phospholipid headgroups, the 
two side minima to the water layer, and the central minimum 
to the terminal methyl moiety of the hydrocarbon chains. 

The bilayer d-spacing and thickness d,, are shown in Fig. 
5 as a function of P/L and at different hydration pressures. 
For all P/L samples at the two different R.H. values in the 
fluid state, d was obtained from the reflectivity curves by 


TABLE 1 
iodinated and noniodinated SCoV E protein 
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fitting the g, values at each Bragg peak as a function of 
Bragg order (i.e., m) to a straight line. The presence of 
peptides in the lipid induces a shift in Bragg peaks toward 
smaller qg, values as compared to the pure lipid. Labeling the 
peptide with iodine did not change the d-spacing signifi- 
cantly. At fixed P/L ratio, the lamellar repeat distance 
decreases with increasing osmotic pressure, reflecting the 
decrease in water layer thickness. The solid lines correspond 
to third-order polynomial fits, which are simply guides to the 
eye and do not represent any theoretical model. At 98% 
R.H., the d-spacing increases as a function of P/L, whereas 
at 75% R.H. the lamellar spacing is approximately constant 
for P/L < 1/10. Above this value a small difference in 
d-spacing between the iodinated and the noniodinated 
samples (i.e., <1 A) is observed. 

The bilayer thickness defined as the distance between the 
two maxima associated with the phosphorus group dpp as a 
function of P/L was determined from p(z). Note that, dpp is 
unaffected by the normalization of the electron-density 
profiles (34). The two different R.H. series in the fluid state 
are presented in Fig. 5 b. An increase in dpp is observed for 
all SARS-CoV E protein concentrations in the bilayer for 
both R.H. series. Only at high P/L and for R.H. = 98% is 
there a significant effect of the iodine label. For the other 


Summary of experimental results for x-ray reflectivity measurements on DMPC mutilamellar bilayer containing 


Hydration condition* Pit 
R.H. = 98% DMPC 1/500 1/100 1/20 1/10 1/7.5 

Integrated intensity* I Non-I I Non-I I Non-I I Non-I I Non-I 
Ip 1 - 1 - 1 1 1 1 1 1 1 

I, 1.39E-1 - 4.76E-2 - 4.7E-2 4.65E-2 4.45E-2 3.0E-2 1.64E-2 5.0E-2 2.6E-2 
lL 2.80E-2 - 1.37E-2 - 0.81E-2 1.15E-2 6.0E-3 5.2E-3 2.15E-3 9.1E-3 3.5E-3 
I; 4.0E-3 - 3.43E-3 - 1.05E-3 2.41E-3 1.0E-3 8.8E-4 3.7E-4 1.3E-3 4.9E-4 
ly 8.0E-5 - 1.7E-5 - 5.86E-5 16.0E-5 6.9E-5 5.6E-5 2.6E-5 7.6E-5 3.7E-5 
Is 3.30E-5 - 5.4E-5 - 1.18E-5 0 0 0) 0 0 0 

Ig (0) - 0 - 0 0 0 0 0 0 0 
R.H. = 75% 

Ip 1 1 1 1 1 1 1 

I, 3.29E-2 7.A1E-3 5.5E-3 2.4E-3 4.2E-3 9.5E-4 2.9E-3 
lL, 1.89E-2 4.6E-3 3.5E-3 2.0E-3 1.8E-3 1.4E-3 1.7E-3 
I; 1.8E-2 4.2E-3 3.6E-3 1.5E-3 1.6E-3 1.3E-3 1.6E-3 
ly 3.35E-5 5.3E-6 3.2E-6 0 0 0 0 

Is 6.98E-4 9.5E-5 4.3E-5 0 2.6E-5 0 3.2E-5 
Ig 3.63E-4 1.7E-6 0 0 0 0 0 

RH. = 11% 

Ip 1 1 1 1 1 1 1 

I, 1.43E-4 4.9E-3 4.9E-3 5.1E-4 3.4E-3 3.9E-4 2.4E-3 
lL, 17.64E-5 0 1.9E-3 0 3.1E-3 0 2.7E-3 
Il 4.78E-2 5.2E-3 8.3E-3 2.9E-3 1.3E-3 2.9E-3 5.1E-4 
ly 3.20E-5 0 0 0 0 0 0 

I; 4.48E-4 2.4E-4 4.3E-5 6.5E-5 6.3E-5 1.3E-4 2.0E-5 
Ig 6.7E-6 - - - - 0 (0) 0 0 0 0 


P/L indicates the protein/lipid ratio, and is given for iodinated (I) and noniodinated (Non-I) samples. 


*Samples hydrated through vapor by equilibration with saturated salt solutions indicated by R.H. value. 


The integrated intensities under each Bragg peak of the reflectivity curve. All the peaks were normalized with respect to the first Bragg peak in the 
reflectivity curve. See text for further explanation. 
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samples, the labeling does not change the structural param- 
eters. For each hydration value, the two data sets (i.e., labeled 
and unlabeled) were fitted to a third-order polynomial, which 
again serves as a guide to the eye. With increasing osmotic 
pressure the bilayer becomes slightly thicker. This phenom- 
enon is well known for pure lipid bilayers and corresponds to 
bilayer thickening upon dehydration. The effect of bilayer 
thickening with P/L is in striking contrast to the bilayer 
thinning observed for many a-helical amphiphilic peptides, 
which have created some interest recently due to antibiotic 
activity. For this class of peptides, bilayer thinning drives 
the transmembrane insertion. Above a critical concentration 
(P/L)*, there is a transition from a parallel to a perpendicular 
(transmembrane) conformation (35-37). 

Finally, the electron-density curves of labeled and unla- 
beled protein are compared for constant P/L and R.H. to 
determine the position of the labeled group in the lipid 
bilayer (Fig. 6). The result supports the helical hairpin 
model, as put forward previously on the basis of a limited 
x-ray data set in combination with FTIR results (20). Fig. 6 
shows the electron-density profile at P/L = 1/10 and 1/7.5 
for the labeled and unlabeled peptide at R.H = 98% and 75%. 
The labeled and unlabeled curves were compared as follows. 
The Bragg peak intensities of both curves were normalized 
to the first Bragg peak of the unlabeled sample. The profiles 
were then computed by Fourier synthesis, the zero-density 
value corresponding to the mean density. Finally, the two 
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FIGURE 4 _ The electron density pro- 
files of the pure DMPC bilayer and 
different P/L ratios of unlabeled (a) 
and labeled (b) SARS-CoV E protein. 
The curves were computed from the 
integrated peak intensities taken from 
Fig. 3 (1e., RA. = 98%) by the FS 
method. The normalization procedure is 
described in the text. 


DMPC 
P/L=1/20 
P/L=1/10 
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curves were multiplied by the same scalar so that the 
maximum in the headgroup region of the iodinated electron- 
density curve is 1. 

At high mol % of the iodine, the electron-density dif- 
ference is large enough to locate the position of the iodine (or 
the phenylalanine) with respect to the bilayer center. Exam- 
ining the profiles in Fig. 6, a and b, reveals a small increase 
of electron density in the headgroup region of the bilayer in 
the presence of the iodine label. In the central region of the 
lipid-peptide profile, corresponding to the bilayer hydropho- 
bic core, the changes caused by iodine are small and below 
the experimental accuracy. Subtracting the electron-density 
profiles of lipid bilayers containing iodinated and unlabeled 
SARS-CoV E proteins should result in an effective iodine- 
density profile. Maxima in these profiles can then be used to 
determine the iodine position (z,) from the bilayer center. At 
higher hydration the pronounced maxima in these curves 
indicate the most probable position of the iodine at z; = 16.5 
and 16.6 A for P/L = 1/10 and 1/7.5, respectively. This 
implies that the iodine-labeled phenylalanine is located ad- 
jacent to the headgroup region. The iodine electronic dis- 
tribution at 75% R.H. (see Fig. 6, c and d) shows a different 
behavior. The oscillatory profile with its small amplitude 
suggests a more disordered state. A small increase in its 
density profile outside the headgroup region could suggest 
that the protein partly aggregates between the bilayers, but it 
is unclear whether this small maximum is really significant. 


I FIGURE 5 (a) The membrane repeat 
distance d as a function of P/L for 
SARS-CoV E protein at T = 45°C and 


. ; per ies i : two different R.H. (b) The bilayer thick- 
fit to 75% data ness dpp, defined as the peak-to-peak 
a mL. R.H=98% distance in the electron-density profiles. 
» |, RH.=98% The iodinated samples are indicated by 
fit to 98% data 
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Anomalous reflectivity 


The above results were solidified by performing an anom- 
alous x-ray reflectivity study (AXR) on the same samples. 
Here, we briefly describe the experimental procedure and 
give a comparison between both sets of results. The ex- 
perimental details were given in a previous article (38). The 
experiment was performed at the European Synchrotron 
Radiation Facility/ID1 (Grenoble, France) by choosing the 
incoming x-ray energy in the vicinity of the Ly, adsorption 
edge (i.e., Ey, = 4.5545 keV) of the iodine. The anomalous 
effect at the iodine adsorption edge is then used to determine 
the label position within the lipid bilayer. Two samples at 
T = 45°C and R.H. = 98% were used in the experiment: pure 
DMPC and P/L = 1/10 with iodinated Phe-23. Reflectivity 
curves were measured at five different energies (4.3975 keV 
< E < 4.5675 keV) on both samples. AXR was first per- 
formed on the pure DMPC sample under the same instru- 
mental and environmental conditions as a test experiment 
and to check for x-ray radiation damage. In fact, no severe 
damage to the sample was observed despite the relatively 
high absorption coefficient at these photon energies. Next, 
the sample with protein was measured at the same five ener- 
gies. As before, all reflectivity and rocking curves indicate a 
highly organized multilamellar film on solid supports. For 
DMPC, more than seven lamellar reflectivity peaks of re- 
flection were obtained, whereas only four to five orders of 
reflection were observed in the presence of SARS-CoV E 
protein due to lamellar disorder and fluctuations (data not 
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FIGURE 6 _ (a and b) Electron-den- 
sity profiles of the P/L: 1/10 and 1/7.5, 
respectively, of the iodinated and non- 
iodinated SARS-CoV E protein at 
R.H. = 98% in arbitrary units. The 
electron-density profiles were normal- 
ized with respect to the first integrated 
peak of the labeled curve. The iodine- 
z(A) density profile (the lower curve in the 
panel) was calculated by subtracting the 
electron densities of the noniodinated 
from the curve of the iodinated SARS- 
CoV E protein. The difference curve 
depicted in the figure has been multi- 
plied by a factor of 3, and the arrows 
indicate the position of the iodine label 
adjacent to the headgroups. (c and d) 
The same profiles as in a and b but at 
R.H. = 75%. The position of iodine is 
not well defined in this case. 
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shown). After plotting the curves against qg, the differences 
can be attributed exclusively to the anomalous effect, as- 
suming that the sample state remains the same over the 
course of the experiment (several hours). Next, all peaks of 
the reflectivity curves are analyzed by the FS method to 
obtain the electron-density profiles p(z, E) for all energies as 
described in Materials and Methods. The iodine-density 
profile was calculated from the differences in the electron- 
density profiles, e.g., Ap;(z, EZ) = pg(z) — pg, (z) (Fig. 7). In 
line with the previous results, the maximum in the iodine 
difference curves is found to be in the headgroup region at 
18.1 A from the center of the lipid bilayer (Fig. 7, arrows). 
However, the density difference curves taken at different 
photon energies agree only qualitatively and do not follow 
the quantitative scaling of the atomic form factor, as has been 
tentatively ascribed to systematic errors in the experiment 
(in particular sample drift in temperature or humidity) (38). 


Effect of peptide concentration on chain ordering 


The effect of adding SARS-CoV E protein to lipid bilayers 
leads to changes not only in the bilayer electron-density 
profile, but also in the short-range ordering of acyl chains. 
Although in most cases adsorbed or inserted peptides and 
proteins further reduce the correlations in the fluid phase, 
leading to a decrease and broadening of the chain correlation 
peak, SARS-CoV seems to have the opposite effect (Fig. 8). 
The scattering distribution as measured in grazing incidence 
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diffraction geometry is shown in Fig. 8 a for both DMPC and 
SARS-CoV E protein at P/L = 1/10 (iodinated) as a function 
of lateral momentum transfer qj). Similar curves were ob- 
tained in the case of the noniodinated protein, and for the 
entire P/L sample series. With increasing P/L, the fluid 
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correlation peak sharpens and moves to higher q). At the 
same time, another component of lipids seems to produce a 
weak maximum at smaller qj), which is less pronounced. 
Note that wide-angle scattering of thin lipid films at in-house 
sources is difficult, yet the general trend is already very clear 
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Same as in panel a but at high hydration 
and temperature (T = 45°C and R.H. = 
98%). A pure DMPC curve is compared 
to curves of different lipid/protein mix- 
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from the data presented. The scans were taken using a home- 
built diffractometer with a sealed tube of Cu Ky (ie., A = 
1.54 A) radiation, equipped with a collimating x-ray multi- 
layer, motorized slits, and a fast scintillation counter. The 
setup is described in detail elsewhere (33). The samples were 
placed horizontally at the bottom of the chamber, and the 
temperature was set to 28°C. 

Lorentzian fits to the scattering distribution data were 
performed for each P/L (i.e., see solid line in Fig. 8 a) to 
quantify the peak position go(P/L) (not shown), correlation 
length €,(P/L) = 1/HWHM (half width at half maximum), 
and, consequently, the nearest-neighbor distance of acyl 
chains dy (i.e., 27/qo). A graph of these quantities as a 
function of P/L is given in Fig. 8 b. It becomes clear that 
SARS-CoV E protein affects the ordering of the acyl tails in 
a dramatic way, even at relatively small concentrations. This 
implies that the protein changes the state of the bilayer in 
its vicinity over some range, and not only locally. Upon 
increase of the protein concentration, the peak shifts to larger 
o Values (e.g., smaller average next-neighbor distance d,y), 
indicating a better ordered packing of the acyl chains. How- 
ever, the observed shift in peak position and the increase in 
€,(P/L) shown in Fig. 8, a and b, could also be due to a shift 
of the main phase transition temperature. It is not clear that 
all samples, in particular those at high P/L, are in the fluid 
state at these temperatures. Since the chain correlation 
maximum at higher temperature deep in the fluid phase was 
hardly visible at the in-house diffractometer, we performed 
additional experiments using synchrotron radiation, for in- 
tensity reasons and better signal-to-noise ratios. These ex- 
periments were carried out at the D4 station of the Doris III 
storage ring of HASYLAB/DESY Hamburg, using a photon 
energy of 20 keV. Importantly, in contrast to Fig. 8, a and b, 
the samples were now heated up in the chamber to the same 
temperatures as in the reflectivity measurements to really 
ensure that the samples were all in the fluid state. In addition, 
the higher water-swelling states were achieved (periodicity 
d ~ 59A for pure DMPC). A huge decrease in the scattering 
intensity with P/L is observed even at small protein con- 
centrations, pointing toward a strong disordering effect of the 
SARS-CoV E protein on acyl chain ordering in the fluid 
phase, see Fig. 8 c. For comparison, the measurements were 
tested at two levels of hydration, one imposed by a pure 
water reservoir in the chamber, and one by a saturated NaCl 
solution. Fig. 8 e shows an approximately linear decrease of 
the nearest-neighbor distance of acyl chains as a function of 
P/L, which is stronger at higher hydration, but much smaller 
than in Fig. 8 b. This comparison supports the assumption 
that the higher P/L samples at T = 28°C in Fig. 8 b were 
already in the gel phase. Note that it is not uncommon that 
the main phase transition changes with P/L. It is also im- 
portant to point out that the changes in peak position are 
small (in the fluid phase) compared to the peak width. 
Finally, Fig. 8 shows the corresponding correlation lengths 
computed from the half width at half maximum values. The 
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samples hydrated from salt solution show an increase in the 
ordering, analogous to the results for T = 28°C in the vicinity 
of the phase transition. Contrarily, the correlation length of 
the highly hydrated samples decreases with P/L. In sum- 
mary, the gel and fluid phases behave quite differently. The 
main effect is always a strong decrease in scattering 
intensity. For high temperatures and high hydration a slight 
decrease in the correlation length and the average next- 
neighbor distance is observed. This is the most relevant 
regime, since the SARS-CoV E protein is in the inserted- 
hairpin state under these conditions. 


SUMMARY AND CONCLUSIONS 


The results show that incorporating the SARS-CoV E protein 
in the DMPC bilayers in the fluid phase induces significant 
changes in the bilayer structure. These changes are nonlocal 
in the sense that at small P/L the effect grows stronger in 
proportion to the lipids that locally surround the protein. This 
can be well understood since a local perturbation in the 
bilayer is generally accompanied by a nonlocal strain fields 
that relaxes the perturbation over some length. Such an effect 
has been observed experimentally for the antibiotic peptide 
Magainin 2 (39), and has been predicted theoretically based 
on bilayer elasticity models (40). In general terms, a hy- 
drophobic mismatch of a transmembrane helix can be the 
source of such a perturbation. In the case described here, the 
experimental evidence gives a consistent picture in that both 
bilayer thickening and a decrease in the acyl chain distance 
are observed. These two phenomena occur simultaneously, 
since the free energy associated with the interaction between 
the chains is minimized if the density in the acyl chains is 
preserved. Thus, to increase bilayer thickness, the chains 
must necessarily approach. Near the gel phase they also be- 
come stretched at high P/L, most likely by an increase in the 
main phase transition. 

In light of these results, let us consider possible molecular 
conformations. From previous FTIR results (20), we know 
that the protein helices are oriented perpendicular. From the 
density profiles of the iodinated protein, we know that the 
Phe-23 group is located at the hydrophilic/hydrophobic 
interface of the bilayer, at least at high hydration. This has 
led us to conclude that the protein forms a small alpha-helical 
hairpin. This is supported by the fact that, corresponding to 
26 amino acids, a total transmembrane length of ~39 A 
would lead to an impossibly high hydrophobic mismatch for 
a transmembrane helix. These values would be applicable, 
assuming that the entire hydrophobic part of SARS-CoV E 
forms an a-helix with each amino acid contributing an axial 
length of 1.5 A (41). In general, protein- or peptide-lipid 
complexes are expected to respond to such an energetically 
unfavorable mismatch situation in a number of ways, de- 
pending on the molecular details of the system. The 
polypeptides can tilt or kink when their TM hydrophobic 
length is too long to match the bilayer, thus reducing their 
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effective length. In the opposite situation, a TM helix could 
adopt a more extended conformation or a nontransmembrane 
orientation (42-44). Alternatively, the lipid bilayer could 
respond to the mismatch situation by ordering or disordering 
their acyl chains or changing the bilayer’s curvature (45,46). 
In our case, it seems that both the protein and the membrane 
responded to this unfavorable situation. Namely, the protein 
adopted a hairpin conformation, with two short helices of 13 
amino acids corresponding to 20 A. At the same time the mem- 
brane changed its thickness due to the presence of SCoV 
E protein. 

For illustration, let us sketch different models in the light 
of these results, to determine which configurations could 
possibly encompass both an inserted hairpin structure and 
the bilayer thickening, and which configurations can be ruled 
out. The first model illustrated in Fig. 9 a shows the protein 
in the hairpin confirmation partly spanning the hydrophobic 
core of the bilayer. This would explain the FTIR and x-ray 
results on the protein conformation, but not why the bilayer 
would tend to thicken. It is also unclear in this case how the 
polar C- and N-termini can be encompassed in the hydro- 
phobic core. In Fig. 9 b, the same conformation, with two 
hairpins forming a homodimmer, is sketched. Our experi- 
ments can so far not distinguish between Fig. 9, a and b, but 
the problems associated with Fig. 9 a obviously also apply 
here. The models shown in Fig. 9, c and d, follow from Fig. 
9, a and b, respectively, if the bilayer thickness is allowed to 
thin to match the hydrophobic length of the hairpin, so that 
the polar ends are facing the other side of the bilayer. 
Obviously the required thinning would be too large and can 
be ruled out, just as the required thickness for one long TM 
helix would be too large (not shown). Finally, the thickening 
observed here directly rules out configurations a—d in Fig. 9. 
Finally, Fig. 9, e and f, shows homodimer configurations, in 
which the proteins associate not laterally but vertically. 
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FIGURE 9 Cartoon representation of different possible models of SARS- 
CoV E protein structure and membrane topology. (a) Monomeric E protein 
inserted into the lipid bilayer membrane. This model suggests that the TMD 
partially spans the lipid bilayer as a hairpin. (b) The same as in a but two 
hairpins forming a homodimer. (c) The same as in a, with the lipid bilayer 
thickness adapted to minimize hydrophobic mismatch. (d) The same as b, 
but again with the lipid bilayer thickness adapted to minimize hydrophobic 
mismatch. (e) Vertical association of two hairpins to homodimers. The 
vertical thickness is adapted to the bilayer thickness by intercalation of the 
two hairpins. (f) Similar to e but with the bilayer thickness increased to 
match the length of two vertically stacked hairpins. 
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Again, simple stacking of two opposing hairpins such as the 
configuration shown in Fig. 9 e is unlikely, since the required 
thickness may be too high. Contrarily, Fig. 9 f, with the two 
opposing hairpins turned by 90° and inserted seems to be 
possible based on simple thickness considerations. This 
putative model assumes that the end termini are packed 
together to suppress the unfavored interactions between them 
and the acyl chains. The model would be in agreement with a 
small positive hydrophobic mismatch, which could explain 
the observed changes in the density profile and the acyl- 
chain scattering. Note that the transition between Fig. 9, e 
and f, is somewhat continuous depending on the degree of 
insertion. In conclusion, the reflectivity results, together with 
the chain correlation measurements, would be in line with 
Fig. 9, e and f, with hydrophobic mismatch considerations 
favoring Fig. 9 f. 

Finally, in this study, we did not find any evidence for a 
destabilization of the bilayer. In fact, nonspecular reflectivity 
(diffuse scattering), as mapped by a two-dimensional CCD 
camera in the geometry of grazing incidence small-angle 
scattering, excludes any dramatic increase of bilayer corru- 
gation or fluctuations with P/L. Thus the bilayers in the 
oriented stacks are found to be quasiplanar with the presence 
of the usual thermal fluctuations (47). A more quantitative 
analysis in the next step could quantify the associated 
changes in the elastic constants, such as the bending rigidity. 
However, an instability can already be ruled out from in- 
spection of the raw data. At the same time SCoV E protein is 
known to play a role in virus budding. Therefore, this phe- 
nomenon must be linked to an asymmetric embedding SCoV 
E protein in the external leaflet of the bilayer that may perturb 
the membrane and also change its permeability. In situ experi- 
ments with vesicles and asymmetric embedding of proteins in 
model systems could help to address this question and may be 
helpful in understanding the effect of the SCoV E protein on 
the formation of ion channels and budding. 
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